Nature of crossover from classical to Ising-like critical behavior 
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We present an accurate numerical determination of the crossover from classical to Ising-like critical 
behavior upon approach of the critical point in three-dimensional systems. The possibility to vary 
the Ginzburg number in our simulations allows us to cover the entire crossover region. We employ 
these results to scrutinize several semi-phenomenological crossover scaling functions that are widely 
used for the analysis of experimental results. In addition we present strong evidence that the 
exponent relations do not hold between effective exponents. 
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It is now well established that the critical behav- 
ior of large classes of systems — including uniaxial fer- 
romagnets, binary alloys, simple fluids, fluid mixtures, 
and polymer blends — belongs to the three-dimensional 
(3D) Ising universality class. However, this behavior is 
only observed asymptotically close to the critical point, 
whereas classical critical behavior is (sometimes) ob- 
served at temperatures farther from the critical temper- 
ature T c , before one enters the noncritical background. 
At intermediate temperatures, a continuous crossover oc- 
curs from one universality class to the other. Due to 
the limited extent of the asymptotic (Ising) regime, the 
thermodynamic behavior in this crossover region is ex- 
tremely relevant from an experimental point of view and 
hence has attracted long-standing attention. Despite this 
attention, several fundamental questions concerning the 
crossover are still open to debate. 

First, we note that the vast majority of all relevant 
studies is limited to the one-phase region only. The be- 
havior of the susceptibility (compressibility) in this re- 
gion is described by means of several approaches, which 
all consist of more or less phenomenological extensions of 
the renormalization-group (RG) description of the criti- 
cal behavior. In Ref. [jD three different descriptions have 
been compared. From a nonlinear treatment of <j> the- 
ory at fixed dimensionality d — 3, Bagnuls and Bervil- 
lier obtained functions describing the full crossover 
from asymptotic critical behavior to classical critical be- 
havior. Belyakov and Kiselev Q carried out a first- 
order e-expansion, which was then phenomenologically 
extended to yield the correct asymptotic behavior in the 
Ising regime. Finally, Ref. Q discusses an extension 
by Chen et al. of the work of Nicoll and Bhattachar- 
jee j5|, based on an RG matching technique. All three 
approaches yield rather similar results and suggest (for 
simple fluids) a smooth monotonic crossover. 

Yet, from the experimental side it has not been so 
easy to confirm this picture or even to make judg- 
ments concerning the quality of the various predictions. 
Widespread attention attracted the findings || on mi- 
cellar solutions, in which exponents were observed that 
did not fit into any of the known universality classes. 
This was analyzed by Fisher M in terms of a crossover 



effect, suggesting that a nonmonotonic variation of the 
susceptibility exponent in the one-phase region is an in- 
trinsic property of the universal crossover scaling func- 
tion. This possibility of nonmonotonicity was then es- 
sentially confirmed in Ref. S , although again an empir- 
ical extension of RG theory had to be invoked as well as 
extremely strong higher-order corrections. All authors 
stressed the need for much more accurate experimental 
results and explicit calculations of the crossover behav- 
ior. Indeed, several experiments have been carried out on 
polymer blends |9|,[l0|], which were analyzed in terms of 
the crossover solution by Belyakov and Kiselev. Polymer 
systems offer the advantage that the Ginzburg number 
ruling the crossover can be varied by varying the length 
of the polymer chains. Nevertheless, it proved difficult 
to cover the full crossover region and results for differ- 
ent polymer blends with widely varying chemical prop- 
erties had to be combined. This in turn led to conjec- 
tures concerning an unexpected pressure dependence of 
the Ginzburg number G, which needed to be fitted to 
each system separately jl0| . A determination of effective 
exponents has not been attempted in these studies. More 
recently, Anisimov and coworkers pd[ | focused on the pos- 
sibility of nonmonotonic variations of effective exponents 
in complex fluids, whereas also for polymer solutions a 
sharp, nonmonotonic crossover of the susceptibility expo- 
nent has been reported []l2| . However, it is unclear how 
the results on micellar solutions fit into this picture and 
also the applicability of the crossover solution by Chen et 
al. to systems which exhibit a nonmonotonic crossover is 
subject to some debate. Nonuniversal crossover (depend- 
ing on parameters in addition to the Ginzburg number) 
is suggested jy], but this behavior may be particular 
to complex fluids due to the occurrence of mesoscopic 
lengths in addition to the correlation length. 

In Ref. fl3]] , we presented a numerical study of the 
crossover to classical critical behavior on either side of the 
critical temperature in two-dimensional systems. This 
revealed a strictly monotonic crossover of the susceptibil- 
ity exponent in the disordered phase (one-phase region) 
and a nonmonotonic variation in the ordered phase. Al- 
though this was an interesting finding in itself, it also 
stressed the importance of a study of 3D systems, since 
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only these systems can be compared with the various 
theoretical crossover descriptions and also an even quali- 
tative difference with the two-dimensional case could not 
be excluded. It is the objective of this paper to present 
the results of a major numerical effort to determine the 
crossover in 3D spin systems. These results allow, for 
the first time, a detailed and rigorous test of the above- 
mentioned crossover functions. 

The crossover is ruled by the parameter t/G, where 
t = (T — T c )/T c is the reduced temperature and G the 
so-called Ginzburg number. Asymptotic critical behav- 
ior occurs for t/G <C 1 and classical critical behavior is 
expected for t/G 3> 1. The additional requirement that 
t must be small implies that only systems with a very 
small G (large interaction range) allow an observation of 
the full crossover. In simple fluids, e.g., the crossover is 
never completed before leaving the critical region. Nu- 
merical calculations offer the advantage that G is known 
precisely and can be made arbitrarily small by increas- 
ing the range of the interactions. Thus, we could vary 
t/G over more than eight orders of magnitude, compared 
to four orders of magnitude in experiments on polymer 
blends f|. The large variation of t/G requires the simu- 
lation of systems with very large coordination numbers. 
Until now, this constituted the main bottleneck for ex- 
plicit calculations. However, the advent of a new Monte 
Carlo algorithm (li) for long-range interactions has now 
allowed us to cover the full crossover region. 

We have simulated classical spin systems, consisting 
of 3D simple cubic lattices with periodic boundary con- 
ditions. The systems are described by the Hamiltonian 
TL/UbT — — Kd{Yi — Tj)siSj, where s — ±1, the sum 
runs over all spin pairs, and the coupling depends on 
the distance |r| between the spins as Kd(r) = ci?~ d for 
|r| < R m and Kd(r) = for r > R m . For any finite R m , 
the critical behavior of this model belongs to the Ising 
universality class, but for R m — > oo it will be classical. 
The consequent singular dependence on R m of all critical 
amplitudes has been derived in Refs. (jl],[l(]]. In order to 
avoid lattice effects we formulate our analysis in terms of 
an effective interaction range R [ [l5| . The R dependence 
of the Ginzburg number is given by G — GoR~ 2d ^ i ~ d \ 
such that the crossover occurs as a function of tR 6 /Go- 
By increasing R we can reach the classical regime while 
still keeping t <C 1 in order to stay within the critical 
region. On the other hand, very small values of t imply a 
strongly diverging correlation length and extremely large 
system sizes are required to avoid finite-size effects. Thus, 
we have constructed all crossover functions by combining 
the results for systems with different Ginzburg numbers, 
such that t had to be varied only within a limited range. 
Systems with linear size up to L = 160 (ca. 4 million 
spins) have been simulated, in which each spin interacts 
with up to 8408 neighbors, corresponding to an effective 
interaction range of 9.8 intermolecular distances. 

We first consider the experimentally most widely stud- 
ied case, viz. T > T c . We have calculated the susceptibil- 
ity from the magnetization density m, x' — L d (m 2 ) /UbT , 



and the scaled susceptibility x = kBT c (R)x' ■ In the 
Ising limit, the latter quantity diverges as Cft~^, with 
7 = 1.237 [0 and Cj = 1.1025 @, whereas mean- 
field theory predicts x — 1/t, i.e., 7mf = 1. The leading 
range dependence of the critical susceptibility amplitude 
is given by R^(i-i)/(i-d) JTJQ ThuS; if onc plots thc 

data as a function of the reduced variable t/G, a data col- 
lapse should be obtained for xGq/R 6 . We will now test 
whether this quantity reproduces the predicted crossover 
behavior and how well it is described by various theoreti- 
cal expressions. First, we consider the phenomenological 
generalization of first-order e-expansion results obtained 
in Q. The main significance of this function lies in its 
rather widespread use for the analysis of crossover effects 
in polymer systems. It is given in the following form, 

t/G=[l + n{xG) eh p- 1)/e 

x{(xG)- 1 + [l + K(xG) e ^r-'/ e } , (1) 

where 8 = 0.508 (25) |l^] is Wegner's correction-to- 
scaling exponent and k w 2.333 a universal constant. 
Asymptotically close to the critical point, x = G 7 ~ 1 (l + 
K 7/0)7 K -7/0£-7 anc j f or t/G » 1 x exhibits the limiting 
behavior 1/t. As the curve reproduces the amplitude of 
the mean-field asymptote, the only remaining adjustable 
parameter is a multiplicative constant in G. This is pre- 
cisely the way G is determined in experimental analy- 
ses. Considering x given by ([!]) as the master curve, 
we set G = [Cf^^(l + K 7/f)-7]i/(7-i) ~ 0.1027 such 
that the 3D Ising asymptote coincides with the curve 
in the limit t/G -C 1. On the other hand, we may 
also calculate the exact Ginzburg number from the ex- 
pression given in Ref. which for our model reduces 
to G = 27/(it 4 R 6 ) w 0.27718/i? 6 . The former value 
Go = 0.1027 leads to a precise reproduction of the Ising 
asymptote and hence has been used in the graph of the 
crossover curve, but the exact value of G may yield a 
better description of the overall crossover behavior. Sec- 
ondly, we consider the crossover function resulting from 
the nonlinear RG treatment of Ref. Q . It is presented in 
terms of a phenomenological function x* {t/9o)/do> which 
represents the RG calculations with a relative error of less 
than 10~ 4 0. The merit of this field-theoretic treatment 
is that it has a more solid foundation than Eq. ([!]) . Also 
X* reproduces the mean-field asymptote and thus we ad- 
just the parameter go such that X*(Vffo)/5o coincides 
with x in the Ising regime. As mentioned before, a third 
crossover function for the susceptibility exists, based on 
an RG matching technique. However, we have not in- 
cluded this function in the present discussion, because of 
its very close resemblance to the solution of Ref. Q] . 

Figure ^ shows our results for the susceptibility in the 
symmetric phase (T > T c ) along with the two theo- 
retical curves. One observes that the simulations ac- 
curately reproduce both the Ising and the mean-field 
asymptote (dotted lines). In between, thc data exhibit a 
gradual crossover, just as predicted by both theoretical 
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FIG. 1. Crossover curve for the magnetic susceptibility 
above T c . The numbers in the key refer to the effective in- 
teraction range R. "BK" and "BB" indicate the crossover 
functions of Refs. and H, respectively. 

curves. However, the vertical scale of the graph covers 
ten decades, which makes it difficult to observe subtle de- 
viations. The logarithmic derivative of the susceptibility, 
= — dlnx/rfln constitutes a much more stringent 
test, see Fig. The excellent collapse of the data for 
all our models with widely differing interaction ranges 
on a master curve suggests that an interpretation of the 
crossover in terms of a universal crossover function (i.e., 
described by a single crossover argument t/G) is appro- 
priate here. In this graph, we have also included a third 
calculation by Seglar and Fisher [|o]|lj , which was shifted 
along the horizontal axis such as to reproduce the initial 
deviations of 7^ from the mean-field value. Apart from 
this somewhat arbitrary shift, the curve resembles the 
other predictions and is given by 

7eff = 1 + (7 " lMY)E[\n(t/G)] , (2) 

with E(]ny) — 1/(1 + y 6 ^ 2 ), where e — 4 — d. Our re- 
sults for 7^ display a smooth, monotonic crossover from 
the Ising value 1.237 to the classical value 1. While this 
agrees with the various theoretical functions shown, it 
clearly contradicts the conjecture of Ref. M, according to 
which a nonmonotonic variation of 7^ in the symmetric 
phase might be a property of the universal crossover scal- 
ing function. The good description of the initial increase 
of 7^ upon approach of the critical point is an encour- 
aging result, since G has been set to its exact value (see 
above) in Fig. || and hence there is no adjustable pa- 
rameter in the solution ([!]). However, we note that the 
figure also reveals a remarkable discrepancy between the 
theoretical calculations and our results. Namely, after 
the initial deviation from 7mf 5 the actual curve proceeds 
with a considerably steeper increase than predicted by 
any of these functions. Following Fisher JtJ, we define 
the gradient V = —d-y^/dlogt. Indeed, T reaches values 
as high as 0.84 per decade, in contrast to the theoretical 
functions for which the maximum of T lies in the range 
0.64-0.74. In practice, this implies that the full crossover 
is completed in one to two decades less than predicted. A 
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FIG. 2. The effective susceptibility exponent 7"^ above T c 
along with three theoretical calculations for this q uantity. 
"SF" refers to the first-order e-expansion of Refs. Pjj7| . 

final interesting aspect of Fig. || is the fact that even for 
R = 1 no significant overshoot of the effective exponent 
above the Ising value can be observed. This is remark- 
able, as it is expected J2TJ that the 3D nearest-neighbor 
Ising model exhibits a negative leading Wegner correc- 
tion. As the form of the curve makes a further increase 
of 7gg for smaller values of t rather unlikely, we conclude 
that the actual effect must be very small. 

Now we proceed to the temperature region below T c . 
We have calculated the susceptibility from the fluctuation 
relation x = L d ({m 2 ) - {\m\) 2 )/k B T . We find that the 
the simulation data faithfully reproduce both the Ising 
and the mean-field asymptote and that the crossover 
curve never deviates far from these asymptotes. The 
overall graph is very similar to that for T > T c . To 
our knowledge it constitutes one of the first determina- 
tions of the full crossover function in the ordered phase. 
Indeed, experimentally it is very difficult to measure the 
coexistence curve and hence the crossover function. Also 
from the theoretical side very few results exist. Bagnuls 
and coworkers |22]| have carried out a calculation using 
massive field theory. Just as for T > T c (see above) they 
present their results in the form of an approximative con- 
tinuous function. However, this function is only valid for 
relatively small values of t/G and does not cover the en- 
tire crossover region. Indeed, for t/G large it approaches 
an asymptote with slope 0.97 instead of the classical 
value 1. Thus, it is not possible to accurately describe 
our results with this function. The exponent 7^ follow- 
ing from our data is shown in Fig. ^ and displays several 
noteworthy features. First, also in the ordered phase 
no nonmonotonicity can be observed within the statisti- 
cal accuracy, in contrast with the two-dimensional case. 
Secondly, the increase of the effective exponent upon ap- 
proach of the critical point is even faster than for T > T c , 
with r as high as 0.11. While 7^ gradually starts to de- 
viate from 7mf when t is decreased, the increase of 7^ 
is rather abrupt. We view this effect as the precursor of 
the "underswing" observed for d = 2 jl^] . In the absence 
of further calculations of 7^, we have attempted to de- 
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FIG. 3. The effective susceptibility exponent y~ s below T c 
"App" denotes the approximation discussed in the text. 
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FIG. 4. The quantity 2 — (7 cff + 2/3 e ff) as a function of t/G. 



scribe the data by a phenomenological generalization of 
the exponent crossover function (|2|). We found that the 
expression £7 (In y) = 1/(1 +y) captures the actual behav- 
ior reasonably well (see Fig. |3|) , although the parameter y 
in the denominator might carry an even larger exponent. 

Finally, we have also considered the crossover behavior 
of the order parameter (|m|) for T < T c . The correspond- 
ing exponent /3 e ff turns out to vary monotonically from 
its Ising value 0.3267 to the classical value 0.5. In Fig. |] 
we display an interesting consequence of this behavior. 
Namely, we have plotted the quantity 2 — (7^ + 2/3 c ff), 
which should be equal to the effective specific-heat expo- 
nent a~g if the standard scaling relations hold between ef- 
fective exponents. However, this quantity varies between 
the classical and the Ising value in a strongly nonmono- 
tonic way, which is very unlikely in view of the smooth 
behavior of 7^ and /3 e ff • Thus, we consider this as strong 
evidence for the violation of scaling relations between ef- 
fective critical exponents. Interestingly, such a break- 
down has been inferred from e-expansions over 18 years 
ago |||, but could experimentally not be confirmed. For 
the first time it has now been explicitly demonstrated. 
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